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Abstract 

We study a feedback stabilization problem for a first-order hyperbolic partial differ¬ 
ential equation. The problem is inspired by the stabilization of equilibrium age profiles 
for an age-structured chemostat, using the dilution rate as the control. Two distin¬ 
guishing features of the problem are that (a) the PDE has a multiplicative (instead of 
an additive) input and (b) the state is fed back to the inlet boundary. We provide a 
sampled-data feedback that ensures stabilization under arbitrarily sparse sampling and 
that satisfies input constraints. Our chemostat feedback does not require measurement 
of the age profile, nor does it require exact knowledge of the model. 

Key Words: bioreactor, hyperbolic partial differential equation, sampled control, 
stabilization. 


1 Introduction 

Age-structured models have been used in mathematical biology and mathematical demog¬ 
raphy for a long time. Models of age-structured populations are either discrete (resulting 
in Leslie matrix models) or continuous (which produce the McKendrick-von Foerster equa¬ 
tion); see [HIS]. Age-structured models have also been used in mathematical economics and 
environmental engineering [3]. The study of continuous age-structured models has focused 
on two different research directions, namely, the study of the dynamics of age-structured 
models, and optimal control problems. Optimal control problems for age-structured models 
have been used in [3l El] and optimality conditions were derived in [21 IE] • The ergodicity 
problem has been studied in many works (such as [101 E]) and many results have been 
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presented for one or mnltiple continnons age-strnctnred population models; see for instance 

m- 

Since the McKendrick-von Foerster equation (which is also called the Lotka-von Foerster 
equation) is a hrst-order hyperbolic partial differential equation (PDF), it is reasonable to 
expect recent results on controlling hyperbolic PDFs (such as [H |6l [71 [131 [m [23]) to apply 
to age-structured models. However, much of this literature uses boundary controls. For 
chemostat models (where the dilution rate is the most commonly used control [HI [13 ITS])- 
it is natural to consider stabilization problems where the dilution rate is used to stabilize 
a specihc age prohle. The dilution rate control enters into the PDF directly, not at the 
boundary. Chemostats form the foundation of much current research in bioengineering, 
ecology, and population biology, and are important in biotechnological processes such as 
waste water treatment plants [13 ED] • Stability properties for the dynamics of chemostats 
have been well studied; see [I71[IH] and the references therein. 

In this paper, we study a simplihed age-structured chemostat model without an equation 
for the substrate concentration, i.e., we consider the substrate concentration to be constant. 
This is justihed in the important case where the inlet concentration of the substrate is used 
to control the substrate concentration (or the substrate concentration is slowly varying). 
Two distinguishing features of the control problem we consider are (a) that the PDF has 
a multiplicative (instead of an additive) input and (b) that the state is fed back to the 
inlet boundary; see the system dynamics 0-0 below. Moreover, we do not require that 
we know the age prohle of the microorganism. Only the value of a linear functional of the 
state prohle (i.e., an output) is known at certain times (namely, the sampling times); see 
([^ below. The controller determines the value of the dilution rate, and its key feature is 
the use of the natural logarithm of the output value. See ([^ for our control formula. It 
is a sampled-data feedback for stabilizing an arbitrary positive valued age prohle, i.e., only 
sampled measurements are required, and the control is constant between sampling times. 
The feedback is valued in a pre-specihed bounded interval, to incorporate input constraints. 

Other key novel features of our work relative to the existing control literature for hy¬ 
perbolic PDFs are that we achieve global exponential stabilization for all positive valued 
initial age distributions, with arbitrarily sparse sampling, and that we do not require exact 
model knowledge. Our key stability estimate is in terms of the sup norm of the logarithmic 
deviation of the state prohle from the equilibrium age prohle; see 0. The proof of our main 
result uses the strong ergodic theorem and the connection between hyperbolic PDFs and 
integral delay equations (IDFs) from our prior work [13]. To our knowledge, this is the hrst 
time that the ergodic theorem has been used to solve a control problem. Our simulations 
show good performance of our control under three operating conditions, and so support our 
view that our work would be useful for industrial applications. 

Definitions and Notation 

We use the following notation. Let X C M be any interval and D C M"- be any set. Let 
C^(X;fl) be the class of all continuous functions / ; X —)■ D, and C^iX^Vt) for any integer 
k > Ihe the class of all functions in C°(X; D) all of whose partial derivatives up through order 
k exist and are continuous on X. Let X°°(X;D) be the equivalence classes of all essentially 
bounded Lebesgue measurable functions / ; X — )■ D with norm ||/||oo = ess sup^gj |/(a)|. Let 
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L^(X; f2) be the equivalence classes of measurable functions f : I fl for which ||/||i < oo, 
where ||/||i = Jj-\f{t)\dt. For each a; G M, let [x] be the integer part of x, i.e., the largest 
integer p such that p < x. We let /Coo denote the set of all strictly increasing unbounded 
continuous functions n : [0, cxd) —)■ [0, cxd) such that k(0) = 0. 

For any subset S' C M and any ^4 > 0, we let PC'^([0, A]; S') denote the class of all 
continuous functions z : [0, A] —)• S' for which there exists a finite (or empty) subset B subset 
of (0,^4) such that: (i) the derivative (d 2 :/da)(a) exists at every point in (0,^4) \ B and is a 
continuous function on (0, A) \ P, (ii) all right and left limits of (dz/da)(a) when a tends to 
a point in the set B U {0, A\ exist and are finite. Given any subset S' C M, we let PC^(X; S') 
denote the set of all piecewise continuous functions, i.e, the set of all functions u : X —> M 
for which there exists a (possibly empty) set P C X such that: (i) u is continuous on X \ P, 
(ii) the intersection of every bounded subset of X with P is finite (or empty), and (iii) all 
right and left limits of u{t) when t tends to a point (from the right or from the left) in the 
set P exist and are finite. 


2 Main Result for Controlled Age-Structured Model 

2.1 Statement of Problem and Theorem 


We consider the age-structured chemostat model given by 

^(ta) + = -{M + D{t))f{t,a) (1) 

for all {t, a) G (0, oo) x (0, A) and 

f(t,0) = f k{a)f{t,a)da for all f > 0 , (2) 

J 0 

where A > 0 is any constant, p : [0, A\ —)■ [0, oo) and k : [0, [0, oo) are continuous 

functions, and we assume that k(a)da > 0. The system (|^-([^ is a continuous age- 
structured model of a population in a chemostat. the boundary condition ([^ is the renewal 
condition, which determines the number of newborn individuals f(t, 0) at each time t > 0, 
A > 0 is the maximum reproductive age, p is the mortality function, / is the density of 
the population of age a G [0, A\ at time t > 0, and k is the birth modulus. Given any 
constants Pmin > 0 and Pmax > -Dmin, the variable D G PG° ([0, cxd); [Pmin, -Dmax]) is called 
the dilution rate and is the control. Physically meaningful solutions of ([^-(|^ are those 
satisfying fit, a) > 0 for all (t,a) G [0, oo) x [0,^]. 

We assume that there is a constant D* G (Pmin, such that 
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|fc(a) exp 



p{s)ds 


da , 


(3) 


which is the Lotka-Sharpe equation [1]. Our assumption that there is a constant D* G 
(Pmin, Pmax) Satisfying ([^ is necessary for the existence of a non-zero equilibrium point for 
In fact, for any constant M > 0, any function of the form 


p{s)ds 


f*{a) = M exp I —D*a — 


■a 


(4) 
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for all a G [0,^4] is an equilibrium point for Q-(|^. Therefore, there are a continuum of 
equilibria. This implies that the dynamics Q-Q cannot be made open-loop asymptotically 
stable to an equilibrium with the constant control D{t) = D*, which is another motivation 


for our globally exponentially stabilizing feedback control design; see Section 2.2 for more 
discussion on the equilibria. 

It is natural to try to design a dilution rate controller D for Q that is a function of values 
of the densities of the newborn individuals, i.e., ([^. However, such measurements may not 
be easy to obtain in practice. On the other hand, it is often possible to hnd a continuous 
function p : [0, A] —)■ [0, cx)) such that we can measure 


y{t) = / p{a)f{t,a)da 


(5) 


at each time t > 0. For instance, the case p{a) = 1 corresponds to measuring the concen¬ 
tration of the microorganisms. Given any desired positive constant lower and upper bounds 
Dmin > 0 and T^max > Drain for the controller, any constant T > 0, and any desired reference 
prohle Q for any M > 0, and setting 

y* = [ P{a)f*{a)da, ( 6 ) 

J 0 

we can prove that our age-structured chemostat dynamics 0-0. in closed loop with the 
piecewise dehned control dehned by 


D(t) = max |i:)min,min |z:>max,F** + F Mn || 

for all t G [iT, {i -|- 1)T) and all integers i > 0, satishes a uniform global asymptotic stability 
estimate for all initial functions /□ G PC^ ([0, H]; (0, cxd )) satisfying 

/o(0) = f k{a)foia)da, (8) 

J 0 

which means that we require that /(0,a) = /o(a) for all a G [0,H]. Our main theorem is: 


Theorem 1. Let A, T, Dmm, and Dmax be any positive constants, with T^max > Dmm- Let 
fi : [0, A] —)■ [0, oo) and k : [0, A\ —)■ [0, oo) he any continuous functions, and assume that 
k G PG^([0, A\] [0, oo)) and is not the zero function. Letp : [0, A] —)■ [0, oo) be any continuous 
function such that p(a)da > 0, and D* G (Pmin, D m^A he any constant satisfying the 
Lotka-Sharpe condition (^. Then there exist a constant a > 0 and a function k G /Coo such 
that for each function fo G PC'^([0, A]] (0, oo)) satisfying 0. the unigue solution of 0-0 
in closed loop with 0 with the initial condition fo satisfies 


max 

0<a<A 


\n(f(t,a)/f*(a))\ < exp(—at)K.( max 

' \0<a<A 


ln(/ci(a)//*(a)) 


(9) 


for all t >0. 
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2.2 Discussion on Theorem [T] 

We discuss the structure of the feedback control, as well as several key features that distin¬ 
guish our controller analysis from existing results. 

The tracking error norm |ln {f{t, a)/f*{a)) \ in the statement of the theorem is motivated 
by the fact that 0 -(§ has the restricted state space X = {f E PC^([0, A]; (0, oo)) : 
/(O) = k{a) f {a)da} . In fact, our logarithmic transformation x{t,a) = ln(/(t, a)//*(a)) 
for (t, a) G [0, oo) x [0, A] produces a control system whose state Xt takes values in M and 
has equilibrium a: = 0, which is the usual setting for hyperbolic PDEs. When M = 1, the 
state space in the new variables is 


X G PC'^([0, A]; M) : exp(a;(0)) = / k{a) f*{a)exp{x{a))da 


( 10 ) 


The nonlinear character of the control problem we consider here is illustrated by the nonlin¬ 


earity exp in (10) 


Consider the special case where the function p in our output (|^ is the birth modulus k. 
Then our output is the density (|^ of newborn individuals. Hence, our theorem is a general 
result on output feedback control that includes the special case where the density of newborn 
individuals is being measured. For our inhnite dimensional systems, the state is the function 
ft, i.e., (/t)(a) = fit,a) for all a G [0,H]. Therefore, even if p = k, our output feedback 
is not a state feedback. The motivation for selecting a specihc M > 0 so that the desired 
equilibrium is (|^ is that for the complete chemostat (which also has an equation for the 
substrate), the selected M > 0 maximizes the yield, in the context of anaerobic digestion; 
see [121 Section 2.4]. 

While there is no explicit Lyapunov functional in our proof of Theorem the function 
Voit) = I ln(/(f, 0)//*(0))| acts as a Lyapunov-Razumikhin functional. This can be seen 


by showing that voit) = ln(/(t, 0)//*(0)) satishes a suitable IDE; see (15), 


and (33). 


The function Vq is a Lyapunov-Razumikhin functional for an IDE, instead of a Lyapunov 
functional for a PDE. Our theorem provides bounds for solutions of the PDF in terms of the 
history of Vq, because 


ln(/(t,a)//*(a))| < Ro(^ - a) + / |P(s)-P*|ds 


( 11 ) 


' t—a 


for aW t > A and a G [0,H]; see (84) below. See also [13], which also uses functions of 
the form Vo(f) = | ln(/(f,0)//*(0))| and which then builds a Lyapunov-Krasovskii func¬ 
tional of the form hP(f) = max{exp(ps)Vo(f — s) : s G [0,H]} for a suitable p G M; 
see the proof of [131 Theorem 2.6]. In the original coordinates, this functional is Q{t) = 
max{exp(pa)g(a)| ln(/(f, a)//*(a))| : a G [0,H]} for a suitable function q. 

As noted in the introduction, to our knowledge Theorem is the hrst application of 
an ergodic theorem to solve a feedback stabilization problem. The proof of Theorem has 
several steps. First, we use the ergodic theorem and a transformation to an IDE to obtain 
estimates for ln(p(iT)/p*). The control is designed such that p(f) —)■ y* as t —)■ oo. In the 
next step, we produce an estimate relating lniyiiT)/y*) and ln(/(f, 0)//*(0)), which we use 
in the third step to show that fit, 0) —)■ /*(0) as f —)■ oo. In the next step, we use our IDE 
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transformation to relate ln(/(t, a)//*(a)) for all a G [0,^4] to ln(/(f, 0)//*(0)). Finally, we 
show that fit, a) —)■ /*(a) for all a G [0,^4]. Onr proof reqnires several lemmas, which we 
turn to next. 


3 Background: Uncontrolled Age-Structured Models 


We review the needed background from mm on uncontrolled chemostats, and other 
material from [16], which we use to prove Theorem below. Let A > 0 be any constant, 
let fi : [0, A] —)■ [0, oo) and k : [0, A] —)■ [0, oo) be any continuous functions, and assume that 
kia)da > 0. Consider the initial value problem given by the two equations 


Oz Oz 

= -k'ia)z{t,a) 


for all (f, a) G (0, oo) x (0, A) and 


zit,0) = / kia)zit,a)da 


( 12 ) 


(13) 


for all t > 0, with an initial condition z{0,a) = z^ia) for all a G [0,24]. System (12)-(13) 
is a continuous age-structured model of a population in a closed ecosystem with no control, 
where /i is the mortality function, x is the density of the population of age a G [0, A] at 
time f > 0, and k is the birth modulus. Physically meaningful solutions are those satisfying 
x(f, a) > 0 for all (f, a) G [0, cxd) x [0,24]. The following existence and uniqueness result 
follows from uni Proposition 2.4] and [ini Theorems 1.3-1.4]: 

Lemma 1. For each absolutely continuous function Zq G C°([ 0 , 24 ];M) such that Xo( 0 ) = 
k(a)zo(a)da, there is a unique function z : [0, cx)) x [0,24] —)■ M that satisfies: (a) For 
each t > 0, the function Zt defined by (zt) (a) = zit,a) for a G [0,24] is in ([0,24];M), (b) 
the function <h : [0, oo) —)• ([0,24];M) defined by <F(f) = Zt is continuously differentiable, 

(c) for each t > 0, the function Zt G ([0,24];M) is absolutely continuous and satisfies 
^i(o) = T kia)zt{a)da, and (d) equation (|i^ holds for almost all t > 0 and a G (0,24). 


Moreover, if zo{a) > 0 for all a G [0, 24], then z{t, a) > 0 holds for all [t, a) G [0, oo) x [0, 24]. 

□ 


We refer to x or the function $ from Lemma as the solution of (12)-(13). We also use: 
Lemma 2. If k e PC^ ([0,24]; [0, oo)), then for every zq G PC^ ([0,24]; M) satisfying 

^o(O) = f kia)zoia)da, (14) 

Jo 

the function z : [0, oo) x [0,24] —)■ M from Lemma^ is on 

iS = {(f, a) G (0, oo) X (0,24) : f — a ^ i? U {0,24}} , 
where B is the finite (or empty) set where the derivative of zq is not defined, and it satisfies 


(12) on S and equation 


for all t > 0. Also, 
z{t, a) = exp I — / /r(s)ds ) v{t — a) 


( 15 ) 
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holds for all {t, a) G [0, cxd) x [0, A], where u G ([—^4, cx)); M) f] ((0, cxd); M) solves 

v{t) = / k{a)exp I — / /i(s)ds ) v{t — a)da (16) 


for all t > 0 for the initial condition v{—a) = exp (/“ /i(s)ds) zo{a) for all a G (0, A]. □ 

In the context of Lemmathe function Zt is of class PC^ for every t > 0. The solution of 
(16) is found by differentiating both sides of (16) with respect to t, then applying integration 
by parts on the interval [0, A], and then solving 


dk 

v{t) = k{0)v{t) — k{A)v{t — A) + / —{a)v{t — a)da, 

./n da 


(17) 


where k{a) = fc(a) exp (—/i(s)ds). Recalling that k{a)da > 0, we also dehne the 
continuous functional P : ([0, A]; M) —>■ M by 


= 


/o'^ ^o(q) (la Hs) exp itijl) + D*) dl) ds)da 
// ak{a) exp (- (/x(/) + D*) dl) da 


(18) 


where D* is from the Lotka-Sharpe condition (|^. Recall the following strong ergodicity 
result, which follows from dl Section 3]: 

Lemma 3. There are constants £ > 0 and K > 1 such that for every absolutely continuous 
function zq G ([0, A];]R) that satisfies and (14), the solution of (12)-(13) satisfies 


exp(J(a))| 2 ;(t, a) — exp {D*{t — a) — J{a)) P( 2 ;o)|da < 
Kexp{(D* — e)t) exp(J(a)) |^o(a)| da for all f > 0, 


(19) 


□ 


where J{a) = IoKs)ds. 

Lemma 1^ follows by choosing {S{t)zo){a) = Zt{a) = z{t,a) for all a G [0,R] as the 
semigroup in m- lik G PC^ ([0, A]; [0, oo)), then for every zq G PC^ ([0, A]; M) that 
satishes 0 , we dehne 


(fit) = exp i—D*t) u(t) — P(^o) for all t > —A, 
where u G (7° ([—^, C)o); M) is the solution of (16) with the initial condition 

u(—a) = exp I / /i(s)ds ) Zoia) 


( 20 ) 


for all a G (OjRj. Then 0 G C°((—A, 0); M). Also, (15) and (19) give 


exp i—D*a) \(pit — a)|da < Pexp (—st) 


exp 


for alH > 0. Therefore, by setting 


C = max^kia)exp ( — / /i(s)ds ) , 


/i(s)ds ) | 2 ;o(a)lda (21) 


(22) 











it follows from (16) and ([^ that 


10(f) I = J k{a) exp ^—D*a — J /i(s)ds ) 0(f — a)da 

< C [ exp (—Zl*a) |0(f — a)| da 

Jq ^ ^ 

< KCexp{—et) / exp I / /i(s)ds ) |^o(a)| da 


(23) 


holds for alH > 0. 


4 Key Lemma 

Our proof of Theorem will also use the following key lemma, which follows from our recent 
results in 


Lemma 4. Let G G 0° ([0, A]; [0, cxd)), set L = G{a)da, and let A G (0, A) be any 
constant such that 


/ G(a)da < 1. 

J 0 

If L > 1, then for each Xq G ([—A, 0); M), the solution x G ([—A, cx)); M) of 

[A 

x(t) = / G{a)x(t — a)da 

J 0 

with the initial condition x{a) = x^^a) for a G [— 24 , 0 ) satisfies 


(24) 


(25) 


minjai, < inf x(t + a) 

-A<a<0 


for all t > 0, where 


< sup x{t + a) < max{a 2 , a 2 &^’'"*'^^} 

-A<a<0 


ai = inf_A<a<oa;o(a), ^2 = sup_^<„<o 2^o(a) , 
c = /o'^G'(a)da, b = ^ , 


(26) 


and h = min{A, A — A}. 


(27) 


□ 


For the proof of Lemma see Appendix |A.l A consequence of Lemma is that for 
every zq G PG^ ([0, A];M) satisfying ^o(O) = k(a)zo(a)da and Zo(a) > 0 for all a G [0, A], 
the solution of (12)-(13) satisfies z{t,a) > 0 for all (f, a) G [0,cxd) x [0, A], To see why, note 
that if 


J k{a)exp^—J /i(s)dsj da > 1, (28) 

then we may apply Lemma l^to ( [T^ , by choosing G{a) = k{a)exp{— fi{s)ds). If 

[ k{a) exp (— [ /r(s)ds^ da < 1, (29) 
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then set x{t) = exp{pt)v{t) for all t > —A for any constant p > 0, where v satishes (16). 
Then 


x{t) = / k{a) exp ( pa 


/i(s)ds ) x{t — a)da for alH > 0 


(30) 


and k{a) exp {pa — ia{s)ds) da > 1 when p is large enough, since k{a)da > 0. 
Also, 

min ( exp ( D*a + / p{s)ds ) Zo{a) ) < min {P{zo) + (j){t + a)) 


0<a<yl 


-A<a<0 


< max {P{zo) + Mt + a)) < max ( exp ( Zl*a + / a(s)ds 2;o(a) 

-A<a<0 0<a<A ' ' ! r y , , y , 


(31) 


for all f > 0. The inequalities ( [M| are obtained by using Lemma with L = 1, x{t) 
from Lemma 1^ taken to be exp(—and G{a) taken to be the integrand in curly 
braces in ([^. Moreover, using our choice (|^ of /*, our formula (15) for the solutions z{t,a) 
of the uncontrolled chemostat (12)-(13), our formula (20) for (j){t), (31), and the fact that 
the solution of the controlled chemostat dynamics ([^^) with D{t) = D* and any initial 
condition /(0,a) = /o(a) with /q G PC'^([0, A]; M) satishes 


/(f, a) = exp{—D*t)z{t, a) 


(32) 


for all (f, a) G [0, cx)) x [0, A] where z{t, a) is the solution of (12)-(13) for the initial condition 
z(0, a) = fo{a), we obtain the following inequalities for all f > 0: 


mina6[a,A](/o(a)//*(a)) < min„g[„,A](/(^, a)//*(a)) 

< max„6[„,^](/(f,a)//*(a)) < mma(.[a,A]ifo{a)/f*{a)). 

The preceding inequalities show that every equilibrium prohle (|^ for every choice of the 
constant M > 0 is stable. However, since every neighborhood of an equilibrium prohle (in 
the norm or in the sup norm) contains inhnitely many equilibria, each equilibrium prohle 
is stable but not asymptotically stable (neutral stability). 


5 Proof of Theorem [T] 


Existence and uniqueness of the solution of the closed-loop system ([^-(|^ with the control 
([^ can be established by the method of steps, as follows. First notice that by Lemma the 
solution z{t,a) of (12)-(13) with the control ([^ and the initial condition zq = /o exists for 
all f > 0, and Lemmas and [^guarantee that Zt is of class PC'^([0, A]; (0, oo) for all f > 0. 
Assume that the solution of Q-(|^, in closed loop with ([^, is dehned on [0,iT] for some 
non-negative integer i and that ft G PC'^([0, A]; (0, cx))) for all t G [0,iT]. Then D{t) can be 
dehned uniquely by ([^ on [iT, {i + 1)T), and D is of class PC°([0, {i + 1)T); [Pmin, -Dmax])- 
Moreover, the solution / of ([^-([^ with the control ([^ satishes 


f{t, a) = exp J D{l)dl^ z{t, a) 


(33) 
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for all { t , a) G [0, cxd) x [0, A] wherever the solution / is dehned. Hence, we are in a position to 
uniquely dehne f{t, a) on [iT, (i + 1)T] x [0, A]. Notice that ft is of class PC'^([0, A]; (0, cxd)) 
for all t in [0, {i + 1)T]. We can continue this process to conclude that the solution of Q-([^ 
with the control ([^ is defined for alH > 0 and satisfies ft G PC^([0, A]; (0, cxd)) for all f > 0. 

Using the fact that the solution of ([^-@ with ([^ satisfies (33) for all {t,a) G [0, oo) x 
[0, H], our choice (|^ of the output gives 


y{t) = exp f— f* p(a)z(t, a)da 


= exp 


Jg D(l)dlj f^p(a)exp (— /i(s)ds) v(t — a)da 
= exp ^D*t — Jg P(/)d/j 

X f^p(a) exp (^—D*a — p{s)ds) (P(/o) + 4>{t — a)) da 


(34) 


for all t > 0, by our choices of v and f from (16) and (20), and the relationship (15) between 


and V. Using (|^ and 0 , we conclude that 

y* = p{a)f*{a)da = where 

fd = f(f p(a) exp [—D*a — fi{s)ds) da . 


(35) 


Combining (34) and (35) gives the following for alH > 0 and t G [iT, {i + 1)P): 
In 


f AA 

V y* 


= \n 

y 


+ D*{t- iT) - P(/)d/ + In 


/ P{fo)+ff g{a)(j){t-a)da 
V P{fo)+ff g(a)0(*T-a)da 


(36) 


= \n 

y* 


( D — n* 1 U — 'iT'] -L ]n f 9A)4‘d 

* V'P(/o)+/o'^ g{a)ct>{,iT-a)da 


where 


g{a) = jd ^p(a) exp I —D*a — 

V ^0 

for all a G [0, H] and 

Di = max i Pmin, min i P^ax, D* + T~^ In 


/i(s)ds 

ym 


(37) 


(38) 


, and 


(39) 


for all integers i > 0, and where the second equality in (36) followed from the sampling 
structure of our controller ([^. We now set 

x{t) = In , Xi = \n (^y^, 

_ / P(fo)+ff gia)4>(ii+l)T-a)da ' 

* V P(.fo)+ff g{a)4>{iT-a)da 

for all integers i > 0. Choosing the positive constant 

6= D*)T,{D* - D^t^)T} , (40) 

we use the following claim: 
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Claim 1. The inequality 


holds for all integers z > 0. 


|a;i+i| < \xi\ - niin{|xj| ,25} + \ui 


(41) 

□ 


For the proof of Claim see Appendix A.2 We also require the following two claims, 
which we also prove in the appendices: 

Claim 2. For all integers z > 0, the inequalities 


Xi > min{0,a;o + z(Zl*-Zimin)?'} 


- min (In 

k=0,...,i 


( P(fQ)+f<^ g(a)0(ir-a)da 


V-P(/o)+/o g{a)4>{kT-a)da 


and 


Xi < max {0, Xo - z (Umax --D*)?"} 


(42) 


max 

k=0,...,i 


Z Pjfo) /o%(a)5h'r-a)da 

W(/o)+/o'^ g{a)<j>{kT-a)da 


are satished. 

Claim 3. The inequalities 

mm{0,x(0)} + mm,.„.(in ( whZltl 

< x(t) < max {0, x{0)} + max^.,.(in 


□ 


(43) 


and 


\x{t)\ < |x[t/T]| 


In 


Pjfo) + fo - a)da 

.Pifo) + 9 {a)<P{[t/T]T - a)da^ 


hold for alH > 0. 


(44) 


□ 


We can combine estimate (43) with our bounds ( [3l| ) on P{zo) + (j){t + a), our choice (|^ 
of /*, and our choices of g and {3 in (35) and (37) (which imply that g{a)da = 1) to get 


k(t)| < |xo| + In 


f maxo<a<A (/o(a)//*(a)) 


Vmino<a<x4 ifoia)/f*{a)) 


(45) 


for all f > 0. The proof of (45) uses the fact that the upper and lower bounds in (31) are 
independent of t. By (23) and (31), we have 


|0(Z)| < K* exp{—et) j /o(a)da, where K* = ZFCexp (^j /z(s)ds^ , 


(46) 


and where K and £ are from Lemma and C was dehned in (22). Let j be the smallest 
integer in [[A/T] + 1, cx)) such that 

A-||/„ll,exp(-a(,r-A))<'|^P(/„), 


where 5 is from (40). We need the following claim, which we prove in Appendix A.5 


( 47 ) 
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Claim 4. For all integers i > j, we have 


Ui 


< 6 and \uA < 


K* II/o111 (exp(^) + 1) exp(£A) 

P{fo) 


exp(—eiT). 


(48) 


Also, 


/ PjM+jA" g{a)<l>{t-a)Aa \ 
V-P(/o)+/o"^ g{a)4>{iT-a)Aa ) 


< 


-^*ll/o|ll(exp(<5)+l)exp(£n) 

Pifo) 


exp {—eiT) 


holds for all integers i > j and all t > iT. 


We next show that 


(49) 

□ 


exp (|a:i+i|) - 1 < exp(-5) (exp (|xi|) - 1) + exp (|nj|) - 1 (50) 


holds for all i > j. When \xi 


< 26, we can use ( |4l| ) to get |a;i+i| < 

\Xi 




which implies 


(50). On the other hand, when \xi\ > 25, we conclude from (41) from Claimj^that |xj+i| < 
\xi\ — 26 + \ui\. The previous inequality, in conjunction with the fact that \ui\ < 6 for all 
i ^ j (which follows from Claim gives 


expd^i+il) - 1 < exp (|a;i| - 25 + |nd) - 1 

= exp (luil) ~ 1 + exp {\ui\ — 26 ) (exp (|a:j|) — 1 + 1 — exp(25)) 

< exp (Iwil) - 1 + exp {\ui\ - 26 ) (exp (|a:i|) - 1) 

< exp (luil) - 1 + exp (-5) (exp (|a;i|) - 1). 


Hence, (50) holds for all i > j. Using (50) and induction, it follows that 


exp(|a;i|)-l < exp (-5(i - j)) (exp (|a;j |) - 1) 

+ (“'*(“ -1-0) (exp(|ti,|) - 1) 


(51) 


holds for all integers i > j. 

Using our upper bounds (48) on \ui\ from Claim|^and the fact that exp(p) — 1 < pexp(p) 
for allp > 0, we get exp (|Mj|) —1 < exp(5) \ui\ for all i > j, and also the following consequence 
of ([5T| for all i > j: 


exp (|xi|) - 1 < exp {-5{i- i)) (exp(|a;j|) - 1) 

_^ X1|/o|5exp(g4)(exp(^)+l) ^ -2-1)) exp(-£/T) . 

1=3 


Since x < exp(x) — 1 < xexp(a;) holds for all x > 0, we conclude that the following holds 
for all i > j: 


Xi 


< exp (^-5(i - j)) expdxjl) |xj| 

^ jCq|/o|5 expM+2,5)(exp(,5)+l) ^ 

1=3 


(52) 
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where 

5 = min{(5, eT}. (53) 

Since S < eT, it follows that exp(—— /)) exp(—e: IT) < exp(—5 i) for all / = j,i — 1 and 
thus (52) implies the following inequality for all i > j'. 


(54) 


|xi| < exp 1^—(5(/— j) j exp (l^jl) l^jl 

^ X1|/o|5expM+2^)(expffl+l) ^^ 


Notice that (54) holds for i = j as well and consequently, (54) holds for all i > j. Using 
(45) and (54) and the fact that x{iT) = Xi for all integers i > 0, we obtain the following 
inequality for all i > j'. 

\xi\ < ^ II/oM^^pW+i) exp(2(5 + eA){i — j) exp ^zj 


+ -p(-U-4)) (W+ln(SSfli» 


(56) 


Since j is the smallest integer in [[^d/T] + l,oo) that satishes (47) it follows that either (i) 
j = [A/T] + 1 or (ii) j > [A/T] + 1 and 


/nl|/„||,exp(-T((j-l)T-yl)) > 


exp(5)+ 1' 


(56) 


In either case, we have 

/i ll/o|li\ exp((5) + 1 

max <1, '■ 


P(/o) j exp(^) - 1 


exp(£(v4 + T)) > exp(ejT) > exp{6j), (57) 


by our choice (53) of S. 


Using (55)-(57) combined with the fact that 

{i — j) exp{—6i/2) < i exp(—5 i/2) < 2exp(—1)/5 


(58) 


for all integers i > j > 0 (which follows because rexp(—r) < exp(—1) for all r > 0), we get 
the following inequality for all i > j: 


\xi\ < max < 1, 


Pifo) 


G\ S{xoJo)\n{S{xoJo)) + 1 exp (-5i/2] , 


where 


G = (exp((5) + l) exp(£yl) max 1, | exp(2(5 - 1), l| and 


S(xq, /o) = ^ 

V KI 1 JUJ mino<a<xl 


maxo<a<A(/o(a)//*(a)) 


l/OW/J /I |\ 

(/o(ci)//"(a))® P'l “IJ' 


(59) 


(60) 


Using (44), (59), (60), the conclusion (49) from Claim|^ and our choice (53) of h, we obtain 
the following inequality for all t > jT: 


|a;(f)| < max < 1, 


K*\\foh 

Pifo) 


G S'(a:o,/o)hi(S'(a:o,/o)) + 2 exp 


(61) 
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Using (45), (57) and (60), we get the following for all t G [0,jT]: 


w*)i < i^oi + In ( rsgti:;:!;;:!:!! j=/»)) 


< exp 


5 [ t' 


[^] ) exp (fj) ln(^(xo,/o)) 


(62) 


< maxjl, I Gexp f-f [h] j ^(xo,/o) In (5(a;o,/o)) 


Estimate (62) shows that inequality (61) holds for all f > 0 


Dehning 


cr = ^ and G = G exp 


(63) 


and using the fact that ^ — 1, we can use (45), (60), and (62) to obtain the following 

for all f > 0: 

|a;(f)| < min jmax |l, } GN'ixo, fo) exp {-2at) , In (^(xo,/o))} , 

where A/'(a;o,/o) = *S'(a;o,/o) In (S'(a;o,/o)) + 2. It now follows directly from the fact that 
min{a, b} < a/^ for all a > 0 and b > 0 that 


\x{t)\< (^max ^GAf{xo, fo) In {S{xo, /o)) j exp (-af) 

for alH > 0. Using ([^ and i§, we get 

y* min (/o(a)//*(a)) < y{0) = [ p{a)fo{a)da < y* max (/o(a)//*(a)) 
0<a<A Jq 0<a<A 


1/2 


(64) 


(65) 


Using dehnition (39) and (65) then gives 


|xo|<lnfmax| max (/o(a)//*(a)),—^- , „ , . ,, , .. 

' V ^ ^^’mmo<a<A(/o(a)//*(a)) 

Next, we dehne the functions 

Qifo) = 

c.<Aifoia)/f4a)) {^aXo<.<A (/o(a)// («)) , mino<.<.(/o(a)//qa)) 


( 66 ) 


mmo<a 


(67) 


and 


«(/„)= (max {1, ) G ( «(/„) In ((3(/„)) + 2 ) In ( (3(/„) 


1/2 


( 68 ) 


It follows from (64) and our formula for S from (60) that the following holds: 


\x{t) I < R{fo) exp (—a t) for all f > 0 


(69) 
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Using our formula for our control -D(t), definition (39), and ( |69| ), we obtain 

\D{t)-D*\ < T-^R{fo)exp{-aT[t/T]) < T-^R{fo) exp {-a{t - T)) (70) 


for all t > 0, because t/T > \t/T] — 1 > {t/T) — 1 for all f > 0. Also, our relationship (15) 
between v{t — a) and the classical solution, combined with our formula (20) for 0(f) and our 
relationship (33) between fit, a) and the solution z{t,a) for the corresponding uncontrolled 
dynamics give 


P{fo) = exp{-D*t)v{t) - 0(f) 

= exp(—i7*f)^(f, 0) — 0(f) = exp{—D*t)exp fit, 0) — 0(f) 

for all f > 0. Hence, our output yit) satishes 


f/(^) == j^p{a)exp{-jlD{t)di^z{t,a)da 
= Cp(a)exp (^- f*D(£)d£ - /^V(s)ds) v(t - a)da 

= /o^P(a)exp J^^B(£)d£ - p(s)dsj {(j){t - a) + P(/o)) exp{D*{t - a))da 

= p{a)exp (^— D{£)d£ — iJ,{s)ds^ (0(f — a) — 0(f)) exp(Zi)*(f — a))da 

+ /o'^P(«)exp (^- Jl; D{£)d£ - p{s)ds^ exp (^-D*a + D{£)d£'^ f{t,0)da 

= f3 f{t, 0) + 0 exp (^D*t - /q D{l)dl^ fo^ ff(a) (</>(£ -a)- 0(f)) da 


for all f > 0, where we used the relationship ( pi3| between z(t, a) and f(t, a), our choice (37) 
of g, the relationship (15) between z(t,a) and v(t — a), our choice (20) of 0, our formula in 
(35) for 0, and the formula (0 for -P(/o)- Dividing through by 0M gives 


/(f,0)/M = 

y(t)/y* - M-^ exp (^D*t - D(l)dl^ 9(a) (0(^ - a) - 0(f)) da. 


Using (46), (70) and the fact that L g(a)da = 1, we get this for all t > A: 


M ^ exp (^D*t — Jg D{l)dl^ f(f 9(a) ((f>(t — a) — (f)(f)) da 

< 2117-^77* exp (-£(f - A)) ll/olli exp (^i?(/o) exp(crT)) 

< 2K* exp {-sit - A)) |ji^ exp {^R{fo) exp{aT)) , 


since 


ll/lloo < M. Using Q, (|33 with /„ = z„, and (|^ gives 


(72) 


(73) 


M ^ exp (^D*t — /J D{l)dl \ g(a) (0(t — a) — 0(t)) da 
<exp(;ifl(/„)exp(<7r)) (74) 

X (maxo<„</i (/o(o)//*(a)) - min (]<„<.4 (Ma)/f’(a))) 


for all f > 0 (by adding and subtracting P(zq) in the integrand in (74)). 
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Combining (73) and (74), we obtain the following for all f > 0: 

M~^ exp (D*t — Jq D{l)<ll\ g{a) — a) — da 


0 Jo 

< exp {^R{fo) exp(crT) - e{t - A)) 


X max ( max (/o(a)/r(a)) - min {fo{a)/f*{a)) ,2K 

I0<a<A 0<a<A J 


l/olli 


(75) 


We next define the following two functions: 


2K* 


V{fo) = exp (^i?(/o)exp(aT)) (^i7(/o) max |i7(/o), 
and H{fo) = maxo<a<A (/o(a)//*(a)) - mino<a<A (/o(a)//*(a)) • 


1/2 


(76) 


Combining (74) and (75), and using the fact that min{a, b} < \/~ab for all a > 0 and 5 > 0, 
gives 

M~^ exp (D*t — f* D(l)dl) g{a) — a) — da 


< V{fo)exp{-e{t - A)/2) 


(77) 


for all f > 0. 


Using (15), our formula (20) for 0(f), our bounds (31) with zq = /o, and our relationship 


(33) between f{t,a) and the solution z{t,a) of the corresponding uncontrolled system, we 
obtain the following for all f > 0: 


(-/o(^(0 - ^min^(/o(a)//*(a)) < 

< exp (^- - D*)dl^ max^ (/o(a)/r(a)) 


im 

M 


(78) 


Using (70) and (78), we obtain the following for all f > 0: 


exp(-^(/o)) min (/o(a)/r(a)) < 

0<a<Jl M 

exp {^R{fo) exp(aT)) (/o(a)//*(a)) 


< 


(79) 


where A{fo) = ^-R(/o) exp(crT). Using dehnition (39) and (69) we get the following for all 


t > 0: 


exp(-R(fo)exp(-trt)) < ^4 < exp (Jt(fo) exp (-cr t)) 


(80) 


Combining (72), (77), (79), (80) and the fact that cr < e/2 (which follows from our choice 
(63) and the fact that b < eT), we obtain the following for all f > 0: 


max<^ exp(-^(/o)) min (/o(a)//*(a)) ,77“(f,/o) 

I 0<a<Jl 

< < min|exp(^(/o))^max^(/o(a)//*(a)) ,77“(f,/o) \ , 


(81) 


where R){t,fo) = exp (i?(/o) exp(-af)) + U(/o) exp (-cr(f-A)). 
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Using the relationship (15) between the function v and the uncontrolled solution z{t,a) 
of (12)-(13) with the initial condition z{0,a) = fo{a) and (33), we obtain v{t) = z{t,0), and 
therefore: 


f{t,a) = exp (^-D{i)di - fi{s)ds^ v{t - a) 

= exp^-D{i)di- p{s)dsj z{t - a,0) (82) 

= exp D{i)de - /i(s)ds) f{t - a, 0)exp D{i)de^ 

when t > a, and 

f{t, a) = exp D{i)di - /i(s)ds) exp /i(s)ds) fo{a - t) (83) 


when t G [0,a), since v{t — a) = exp{J^ */i(s)ds)/o(a — t). Hence, we can use the formula 
(|^ for f* to obtain 


f*{a) 


exp 



{D{l)-D*)dl 


M 


for all (t, a) G [0, cx)) x [0, A] such that t > a, and 


f*{a) 




t) 


for all (t, a) G [0, cx)) x [0, A] such that t < a. We now dehne the functions 


(84) 


(85) 


Bi{t,fo) = max{exp(-C(f,/o)) - S(f,/o), 

exp (-^^(/o)exp(crT)) mino<a<A {fo{a)/f*{a))} , 
B 2 {t, /o) = min {exp (U(f, /o)) + S(f, /o), 

exp {-^R{fo) exp{aT)) maxo<a<A (/o(a)//*(a))} , 
C{t, fo) = R{fo) exp{-a {t - H)), and 
S(U/o) = U(/o)exp(-cr(f - 2H)) . 


Combining (81), (84), (85), and (70) gives the following for alH > 0: 


exp (-HT"V“^ exp(-cr(f - T - A))R{fo)) Bi{t, fo) 

< mino<a<A (/(Ua)//*(a)) < maxo<a<A (/(t, a)//*(a)) 

< exp (HT"exp(-cr(f - T - A))R{fo)) B2{t, fo) 


We now set 


w{t) 


max 

0<a<yl 


|ln(/(Ua)/r(a))|, 


Wi{t) 

W2{t) 


max In (/(f, a)//*(a)) 

0<a<A 

max ln(/*(a)//(f,a)) 

0<a<A 


and 


( 86 ) 


(87) 


( 88 ) 
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for all t >0. Clearly, definition (88) implies that 


w{t) = max |ln(/(t,a)//*(a))| 

0<a<A 


= max {max{ln(/(t,a)//*(a)) ,ln(/*(a)//(t,a))}} 

0<a<A 


= max ( In (/(t, a)/f*{a)) , max In (/*(a)//(t, a)) 
I 0 <a<A 0 <a<yl 

= In (max | max (/(t, a)/f*{a)), max (/*(a)//(t, a)) 

\ I 0<a<n ^ A 


(89) 


0<a<A 


for alH > 0, from which we get 


max (/(t,a)//*(a)) < exp(w(t)) 

0<a<A 

and min {f{t,a)/f*{a))>exp{—w{t)) 
0<a<A 


(90) 


for all f > 0. Therefore, our dehnitions (18), (67), and (68) for P, Q, and R give 


P{fo) > exp{-w{0))P{f*), Q{fo) < exp(3u;(0)), R{fo) < Q{w{0)), and 
ll/olli < exp(M;(0)) ||/*||^, 


where 


Q{s) = 3exp(3s) 


max 


K*\\f% 

P{f*) 


1/2 

((s+l)®)'-'" 


(92) 


for all s > 0. 


Also, our definition of V{fo) in (76) in conjunction with our bounds (90) and (91) give 


!/(/„) < P(.w(0)), 


(93) 


where 


p{‘>) = 

exp r^(5(s) exp((TT)') (^2 sinh(s) max |2 sinh(s), gxp(s(; 


ll/‘ll 


1/2 


(94) 


for all s > 0. Also, our dehnitions of Wi and W 2 in (88) in conjunction with estimate (87) 


and our bound on R{fo) in (91) give the following for all f > 0: 


W 2 {t) < AT ^ exp(—(j(f — T — A))Q(t(;(0)) + In (5]^ ^(f,/o)) and 
wi{t) < exp(-a(f - T - A))Q(w(0)) + In {B 2 {t, fo)) 


(95) 


Dehnitions (86) in conjunction with (90), our bounds on R{fo) and V{fo) in (91) and (93) 
and the facts that ln(a + b) < ln(a) + a~^b for all a > 0 and 6 > 0 and min{a, b} < \/ab for 
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all a > 0 and b > 0, imply that 


ln(52(t,/o)) < 

min |^exp(crT)(5(w(0)) +w(0),ln (exp (C(t,/o)) + S(t, 

< min |^exp(crT)(5(w(0)) +w{0),C{tJo) + S(t, /o) exp (-C(t, 

< min I^ exp(crT)Q(w;(0)) + w;(0), exp (-cr(t - 2^)) (i?(/o) + V (fo)) 

< exp (—cr(t — 2A)/2) 

exp(trT)Q(w(0)) + w(0)j (p(ui(0)) + Q(w(0)) 


(96) 


holds for all t > 0. Also, (86) in conjnnction with (90) and onr bonnd (91) on R(fo) give 

In /o)) < -y;Q(w( 0)) exp(aT) + w(0) for alH > 0 (97) 

and 

In fo)) < - In (exp fo)) - fo)) (98) 

for all f > 0 snch that 1 > exp {C(t, fo)) S(t, fo). Using the facts that ln(l + x) < x and 
— 1 < e^x hold for all x > 0 and (98), we obtain the following for all f > 0 that satisfy 
1 > 2exp (C(f,/o))S(f,/o): 


ln(Bp((,/„)) < In ( „p(_c(,j‘))_s(.j„) ) 


= In 1 + 


< 


l-exp(-C(f,/o))+S(t,/o) 

exp(-C(i,/o))-S(i,/o) 

exp((7(t,/o))-l+Sp,/o)exp(C(t,/o)) 

l-S(t,/o)exp(C(i,/o)) 


(99) 


< 2 (exp {C{t, fo)) - 1 + S(t, fo) exp {C{t, fo))) 

< 2(C(f,/o)+ S(f,/o))exp(C'(f,/o)) 


Next note that that by onr dehnitions (86) and onr bonnds on R{fo) and V{fo) in (91) 


and (93), the ineqnality 1 > 2 exp (C(f, fo)) H(f, fo) holds if 

1 > 2exp((5(r(;(0))exp(2(TA))P(ta(0))exp(—cr(f — 2A)), 


which holds if 


0 > ln(2exp(Q(r(;(0))exp(2crA)) + ln(P(t(;(0)) + 1) — cr{t — 2A). 


( 100 ) 


On the other hand, (100) holds if 

t > 2A + a~^ ln(P(t(;(0)) + 1) + cr~^ ln(2exp(Q(w(0))exp(2crA)). 
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Consequently, we conclude from (99), (86), our bound on R{fo) from (91), and (93) that 
lii(Br'(i./o)) < 

2exp ("(5(ta(0)) — cr(f — 2^4)') ('p(tc(0)) + Q(«;(0))') for alH > T(rt;(0)), 


where T(s) = 2^4 + cr ^ ln(P(s) + 1) + a ^ ln(2exp((5(s)exp(2(TA)) for all s > 0. Then (97) 


and (101) give: 


In (Pi /o)) < exp (t - T(m;(0)))) 

X max 11(;(0) + ^Q(ta(0)) exp((jT), 2P(w(0)) + 2Q{w 


( 102 ) 


for all f > 0. Also, (89) anc 


our defi nition s of wi and W 2 in (88) give w{t) = max {tci (t ), W 2 


[m 


for all t > 0. Using ( |95[ ), ( |96[ ) and ( |102[ ) and noting that (a) the functions P and Q in ( |94[ ) 
and 


), respectively, are of class /Coo and (b) the function T is non-decreasing, we con¬ 
clude that there is a function n G /Coo such that w{t) < exp {—at/2) n (ta(0)) for all f > 0. 


Therefore, the theorem follows from our definition of w{t) from (88). 


6 Simulations 


To demonstrate our control designs from Theorem [T| we carried out three simulations. In 
each simulation, we took the horizon A = 2, the constant mortality function /u(a) = jJ. = 0.1, 
P* = 1, and the birth modulus 


k{a) 
9 = 


ag, 

(2 - a)g, 


(M+D*)2 

(l-exp(-(M+n>*)))" 


a G [0,1] 
a G [1,2] ’ 

= 2.718728 


where 


(103) 


The constant g is chosen such that the Lotka-Sharpe condition (|^ holds with P* = 1. The 
output is 

Vit) = Jo f{t,a)da (104) 

which is the total concentration of the microorganism in the chemostat. Our objective is to 
stabilize the equilibrium profile 


f*{a) = exp {—{D* + fi)a), aG[0,2]. (105) 

The equilibrium value of the output is 

y* = f*{a)da = = 0.808361. (106) 

We tested the output feedback law 

^ ('107') 

max{P^in,min{P^a,„P* + T-Mn(/(*T,0)/r(0))}}, t e [tT, (t + 1)T) ^ 
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and the output feedback law 


D(() = 

max min {£>„„,!>•+T Mn (!/(!T)/!/*)}} , t e [iT, [i + i)T) 


(108) 


where i = 0,1, 2,..and for both controllers we chose T = 0.4, Umin = 0.5, and /Imax = 1-5. 
We took our initial conditions to have the form 


fo{a) = bo — bitt + cexp{—9 a), aG[0,2] 
where bo, c, and 6 are positive parameters that we specify below, and where 

bi = g~^{g - l)bo + ce~‘^ (1 - exp(-6'))^ - cg~'^ 


(109) 


( 110 ) 


is chosen so that /o(0) = k{a) fo{a)da holds. We must also choose the parameters such that 
minag[o, 2 ] /o(a) > 0. We generated the simulations using a uniform grid of function values 
f{ih,jh) for i = 0,1,..., 50 and i > 0, where h = 0.04. For i = 0 we had /(O, jh) = fo{jh) 
for j = 1, ...,50, where fo is from (109). 

We computed the integrals 


y{ih) = f{ih,a)da and f{ih,0) = k{a) f {ih, a)da 


( 111 ) 


numerically for each i > 0. Since we wanted the numerical integrator to evaluate the integrals 
(111) exactly for every i > 0 when f{ih,a) = Cexp{aa) for certain real constants C and a, 
we did not use a conventional numerical integration scheme, such as the trapezoid rule or 
Simpson’s rule. The reason we wanted to evaluate the integrals exactly when f{ih,a) is an 
exponential function is that the equilibrium prohle (105) is an exponential function and we 
would like to avoid a steady-state error due to the error induced by the numerical integrator. 
To this end, we set 


>C(q j, h) = In (j + l)h)) - In {f{ih,jh)) 

and I{i) = {j : f{ih, (j + l)h) = f{ih,jh)}, 


and we used the integration schemes 

C{j+l)h 


f{ih,a)da ^ Ii{j) = 


Ijh 


U f{ih,{j+l)h)-f(ih,jh) 
C{i,j,h) ' 

hf{ih,jh), 


] i m 

3 e 1(0 


( 112 ) 


(113) 


for j = 2, 3,..., 49 and ^ > 0, and 


f{ih,a)da ^ /*(!) = 


h 


f(ih,2h)—f^{ih^h)/f{ih^2h) 


ln(f(ih,2h)) — \n{f(ih,h)) 

2hf{ih, h), 


if f{ih, h) ^ f{ih,2h) 
if f{ih, h) = f{ih,2h) 


(114) 


and we set 


af{ih,a)da ^ Ji{j) = 


ijh 


C{i,j,h) 




ifiih,ij+l)h)-f(ih,jh)) 


3 e m 


(115) 
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for j = 2, 3,24 and ^ > 0, and 

af{ih,a)da ^ Ji{l) = 

2h^ f{ih^2h) 


and 


ln(/ (ih,2h))—\n{f(ih,h)) {\n{f{ih^2h))—ln{f{ih,h))) 

2h? f{ih^ 2h), 


( f{ih,2h)— 

iff^,h,h)^f{th,2h) 


if f{ih, h) = f{ih,2h) 


(116) 


j^ 0 +i)/i (2 _ a)f{ih,a)da ^ Ki{j) = 

h ^ hf(ih,{j+i)h)-hf{ih,jh) j ^ 

J e m 


h^mxj+Dh) 


(117) 


C(i,j,h) 


(2 - ^/i) hf{ih,jh), 
for j = 25, 26,..., 49 and i > 0. 


The derivation of formnlas (113)-(117) is based on the interpolation of 


/j(“) = Cjexp(irja) 


(118) 


throngh the points {jh, f{ih, jh)) and {{j + l)h, f{ih, {j + l)h)) for j = 1,2,...,49. More 
specihcally, we obtain the following for j = 1, 2,..., 49: 


aj = h nn{f{ih,{j+ l)h)/f{ihjh)) and 
Cj = f{ih,jh) {f{ih, {] + l)h)/f{ih,jh))~^ 

Using this interpolation, the exact integration formnlas are nsed. For example, for 

r{j+dh 

/ a/(ih, a)da for j = 2, 3,..., 24, 

Jjh 

we get the following when aj = h~^ In {f{ih, (j + l)h)/f{ih,jh)) ^ 0: 


(119) 


( 120 ) 


fjh 


{j+dh 


f{ih,a)da ~ afj{a)da = Cj aexp{aja)da 


= Cjaj hexp(crjjh) ((j + 1) exp(crjh) - j) - Cja^ {exp{aj{j + l)h) - exp{ajjh)) 
On the other hand, when aj = h~^ In {f{ih, (j + l)h)/f{ih,jh)) = 0, we get 

/•{j+dh r{j+dh _ rU+dh. ^2 

/ af{ih, a)da ~ / afj{a)da = Cj / ada = —Cj {2j + 1). 

Jjh Jjh Jjh 2 


Combining the above formnlas with the estimated valnes for Cj and a,- in (119), we obtain 


formnla (115). Similarly, we derive formnlas (113), (114), (116) and (117). Notice that the 


formnlas (113), (114), (115), (116), and ( |117 ) allow the nnmerical evalnation of the integrals 
(111) for every i > 0 withont knowing f{ih, 0). 
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Since the time and space discretization steps are both h, we have the exact formula 


f{{i + l)h,jh) = 

f{ih, {j — l)h) exp {—{fi + Di)h) for j = 1, 2,50 and alH > 0 , 


( 121 ) 


where Di = D{ih). Therefore, we have the following algorithm for simulating the closed-loop 
system: 

Algorithm: Given f{ih,jh), for j = 1,50 and certain i > 0 do the following: 


1. Calculate f{ih, 0) ^ g MJ) + 9 EZ 25 Ki{j), where Ji{j) and Ki{j) are given by 


^49 


(115), (116), and (117). 


2. Calculate y{ih) ~ where Ii{j) is given by (113) and (114). 


3. If Y is an integer, then set 


Di = max {Dmin, min {Umax, 71* +T ^\n{y{ih)/y*)}} ; 
otherwise, set = Di_i. 


4. Calculate f{{i + l)h,jh), for j = 1,...,50 using (121) 


The above algorithm with obvious modifications was also used for the simulation of the 
open-loop system, and for the simulation of the closed-loop system under the output feedback 


law (108). We next present the results of our three simulations. 


In our hrst simulation, we used the parameter values 5o = 0.2, hi = 0.15184212, c = 0.8, 
and 0 = 1 in our initial conditions. In Figure 1, we plot the control values and the newborn 
individual values. We show the values for the open loop feedback D{t) = 1, and for the state 


and output feedbacks from (107) and (108). Our simulation shows the efficacy of our control 
design. 

In our second simulation, we changed the parameter values to 60 = 1, 7i = 0.7592106, 
c = 4, and 6 = 1 and plotted the same values as before, in Figure 2. The responses for the 


output feedback law (107) and the output feedback law (108) are almost identical. 


This second simulation was made with an initial condition which is not close to the equi¬ 
librium profile (in the sense that it is an initial condition with very large initial population). 


The difference in the performance of the feedback controllers (107) and (108) cannot be 
distinguished. 

In our final simulation, we tested the robustness of the controller with respect to errors in 
the choice of D* being used in the controllers. We chose the values bo = 0.2, hi = 0.15184212, 


c = 0.8, and 9 = 1, but instead of (107) and (108), we applied the controllers which are 


defined for all t G [iT, (i -|- 1)T) and for integers f > 0 by 


and 


D{t) = max jZli, 


D{t) = max \^D-a 


. mm 


0.7 + T-1 In U^{^T, 0)/rm} } 


mm 


{£>„„,0.7 + r-‘ln( 9 (jr)/ 9 ')}}. 


( 122 ) 


(123) 
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1.2 n 

D(t) 

1.05 - 


0.9 


0 0.5 1 1.5 2 2.5 t 3 


Figure 1: First simulation. The red line is for the output feedback (107), the blue line is 
for the output feedback (108) and the black line is for the open-loop system with D{t) = 1. 


In both cases, we obtained lim^^+oo fit, 0) = 1.1275 and 

lim /(t, 0) = 1.1275 and lim Dit) = D* = 1. (124) 

t —^~1“00 t —^~|“00 

Hence, a —30% error in D* gave a -1-12.75% steady-state deviation from the desired value of 
the newborn individuals. See Figure 3. 

Notice that a constant error in D* is equivalent to an error in the set point since we have: 


D{t) = max min {Umax, 0.7-FT Mn (/(iT, 0)//*(0))}} 

= max {D^^, min D* + In (/(^T, 0)//*(0)) - T-^0.12}} (125) 

= max {Drain, min {Umax, D* + T-^ In ifiiT, 0)/(1.1275/*(0)))}} 
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0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 f 5 


Figure 2: Second simulation. The red line is for the feedbacks (107) and (108) and the black 
line is for the open-loop system with D{t) = 1. 


for the output feedback case (107) and 


D{t) = max {Umin, min {Umax, 0.7 + T“Mn (j/(iT)/j/*)}} 

= max {Dmin, min (Umax, + T-'^ In {y{iT)/y*) - T“i0.12}} (126) 

= max {Dmin, min (Umax, +T-^ln {y{iT)/{1.1275y*))}} 


for the output 


feedback case (108). An interesting feature of the closed-loop 


lim D(t) = D* = 1. 

t —^~1“00 


system is 


that 

(127) 


It may be worth considering an adaptive strategy for the elimination of errors in D* (i.e., a 
hybrid strategy that adapts the applied value of D*). We leave the search for such a strategy 
for future work. 
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1.25 - 


1.1 - 
D(t) 

0.95 - 

0.8 - 


0.65 


0.5 ^^^^^^-, 

0 0.5 1 1.5 2 2.5 t 3 


Figure 3: Third simulation. The red line is for the output feedback (122) and the blue line 
is for the output feedback (123). 


7 Conclusions 

Chemostats play a vital role in biotechnological applications, such as the production of in¬ 
sulin and in waste water treatment plants. Age-structured chemostats produce challenging 
control problems for hrst-order hyperbolic PDFs that are beyond the scope of the existing 
controller methods for ODEs. We studied the problem of stabilizing an equilibrium age 
prohle in an age-structured chemostat, using the dilution rate as the control. We built a 
sampled-data dilution rate feedback control law that ensures stability under arbitrary phys¬ 
ically meaningful initial conditions and arbitrarily sparse sampling. Our control does not 
require measurement of the whole age prohle, or exact model knowledge. The proposed 
feedback also applies under arbitrary input constraints. The proof of our main result is 
based on (a) the strong ergodic theorem and (b) our approach from [13] for transforming 
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a first-order hyperbolic PDE into an integral delay eqnation. Onr simnlations demonstrate 
the good performance of onr controllers. 

We hope to bnild on onr research, in two ways. First, since the growth of the microorgan¬ 
ism may sometimes depend on the concentration of a snbstrate, it would be useful to solve 
the stabilization problem for an enlarged system that has one PDE for the age distribution, 
coupled with one ODE for the substrate (as proposed in [22], in the context of studying 
limit cycles with constant dilution rates instead of a control). Second, it would be useful 
to extend our work to cases where the control is subject to uncertainties, and then seek 
generalizations of our exponential stability estimate such as input-to-state stability under 
input constraints and sampling. Finally, we hope to cover state constrained problems, which 
add the requirement that the states must stay in prescribed subsets of the state space for all 
nonnegative times, in addition to the nonnegativity requirements on the physical quantities. 


Appendices 

A.l Proof of Lemma |4| 

Local existence and uniqueness for every initial condition Xq G ([—A, 0); M) is guaranteed 
by m Theorem 2.1]. We dehne the functions 


14(t) = sup x(t + a) and W(t) = inf x(t + a) 

-A<a<Q -A<a<0 


(A.l) 


for all t > 0 for which the solution of (|25|) exists. Let q > 0 and t > 0 be sufficiently small 


such that the solution exists on [t,t + q). From equation (25), we get 


V^{t + q) = sup x{t + q + a)= sup x(t + s) 

—A<a<0 q—A<s<q 


= max <( sup x{t + s), sup x{t + s) 

^q—A<s<0 0<s<q 


<max<K(t), sup G{a)x(t + s — a)dc 


0<s<q 


= max 


V^{t), sup <G(a)a;(t-|-s — a)da-|- G(a)a;(t-h s — a)da 

0<s<q 1 


< max < V^(t), sup < sup x(t + 1) G(a)da -|- sup x(t + 1) G(a)da 

0<s<q A</<s—A s—A<l<s 


<max<K(t), sup a;(t-|-/) A G(a)da-|- sup x(t + 1) G{a)da 

1^ —A<l<q—A —A<l<q 

Using the facts that 

J^G{a)da = L>l and /^'^G(a)da<l (A.2) 

and assuming that q < min{A, A — A}, it follows that 

14(tg) < max |l4(t), K(t) — [ G(a)da^ -h 14(tg) [ G(a)dal , (A.3) 



which gives V^{t + q) < max{6i4(^), K(^)} with b as in the statement of the lemma, by 


separately considering the two possible values for the maximum on the right side of (A.3). 


(We need the maximum max{614(t), K(t)} instead of just K(t), to allow the possibility that 
K(t) is nonpositive.) 

Similarly, we can use the decomposition 


W(t + q)= inf x(t + q + a) = inf x(t + s) 

-A<a<0 q-A<s<q 

= min < inf a;(f + s), inf + s) > , 

I q-A<s<Q 0<s<q I 


the definition of hh(t), equation (25), and the fact that q < min{A, A — A} to get 


W{t + q) 

>min<hh(f), inf G{a)x(t + s — a)dc 

I 0<s<q ^ 


= mm 


iy(f),^mf IG(a)a;(f + s — a)da + G(a)a;(f + s — a)da 


>min<W(f), inf < inf a;(f + Z)G(a)da + inf x(t + 1) f G(a)da 

I 0<s<q Is—yl<(<s—A ^ s—A<l<s ^ 

>min|hh(f), ^ inf +/)/^ G(a)da + mf x(t + 1) G{a)dc 


— A<l<q 


> min 


, W{t) (l - G'(a)da) + W{t + q) G'(a)da} , 


so (A.2) gives W{t + q) > mm{bW{t),W{t)}. It follows from (A.3) that the solution of 
(25) is bounded on [t, t + q) when q < min{A, A — A}. A standard contradiction argument 
in conjunction with [T3l Theorem 2.1] implies that the solution exists for all f > 0. Using 


induction and (A.3), we can now show that 


min{&lU(0),lU(0)} < W{ih) < V,{ih) < max{&*K(0), 14(0)} 


(A.4) 


for all integers Z > 0, where h = min{A,A — A}. Inequality (26) now follows from the 


dehnitions of K and W and (A.3) and (A.4), by choosing i such that i < t/h. This proves 
the lemma. 


A.2 Proof of Claim [T] 


We distinguish between the following cases. 


Case 1: /Imin < D* 


+ T ^Xi < 


In this case, our choices of the Hj’s and Xj’s imply 


that Di = D* + T~^x. 

Zlj’s, the fact that Di = D 
we get Xi+i = Ui, which directly implies (|41l) 


Using our expressions (34) for the output y{t), our choices of the 
+ T~^Xi, and the fact (f){t) and y{t) are continuous mappings. 
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Case 2: -Dmin > D* + T ^Xi. In this case, (38) and (39) give Di = 


continuity of 0(t) and y{t) and setting t = 
the fact that = -Dmin, we get 


Using the 

[i + 1)T in (34), and using dehnition (38) and 


Xi+i = Xi + {D* - Drain)T + Ui- 

The inequality Dmin > D* + T~^Xi implies that 0 > (D* — Dmin)?" + Xi and 0 > Xi. Using 
the previous inequalities, the equality Xj+i = x* + (D* — Dmin)T + Ui, and our choice of 5, 
we get 


|a:i+l| < |Xi + (D* - Dmin)T| + |Mi| = -Xi - (D* - Dmin)T + 

< |xi| - 2(5 + \ui\, 


which again gives (41). 


Case 3: D* + T ^Xj > D 


Then (38) and (39) give Dj = Dmax- Since (j) and y are 


continuous, we can set t = {i + 1)T in (34), and use (38) and the fact that Dj = Dmin 
to get Xj+i = Xi — (Dmax — D*)T + Ui- The inequality D* + T~^Xi > Dmax implies that 
Xi — (Dmax — D*)T > 0 and 0 < Xi. Then the equality Xi+i = Xi — (Dmax — D*)T + Ui give 


|Xi+l| < |Xi - (Dmax - D*)T| + |Mi| = Xi - (Dj, 

< |Xi| - 2(5 + \Ui\, 


D*)T+\ui\ 


so 


(41) holds again. This proves Claim 


A.3 Proof of Claim |2l 


The proof of (42) is by induction. First notice that both inequalities in (42) hold for i = 0. 


Next assume that (42) hold for certain integer i > 0. We consider three cases. 

+ T“^Xi < D„ 

Ui 


Case 1: Dmin < D* 
of Claim gives Xi+i = 


In this case, our treatment of Case 1 in our proof 


Consequently, our dehnition (39) gives 


f ^(U)+/o^ g{a)4>{{i+l)T-a)da 
* V P(/o)+/(f g(a)<j!)(iT-a)da 

^ f ^’(/o)+/cf g{a)<l,iii+l)T-a)da 

— maxfc = o ,..., i+i I m ( yj 

V V Pifo)+Jo g{a)(t>ikT-a)da 

< max { 0 ,Xo - (* + 1 ) (Dmax - D*)T} 

, f P{fo)+fi^ g{a)(t>i(i+l)T-a)da 

+ maxfc = o ,..., j+i I m ( — a , , — r ,— 

V V ^(/o)+Jo g(a)4>{kT-a)da 


which implies the second inequality in (42) with i + 1 in place of i > 0. Similarly, we obtain 


the hrst inequality in (42) with i + 1 in place of i > 0. 

Case 2: Dmin > D* +T“^Xj. Arguing as in our treatment of Case 2 in our proof of Claim 
we get Xj+i = Xi + (D* — D m\r, )T + m,. Hence, (39) and (42) and the fact that D* > Dmin 
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give 


^2 + 1 

> 


Xi + {D* — -Dmin) T + Ui 
min {0, Xo ^i{D* - Dram) T} + {D 
+ minfc=o,...,i I In 


mm, 

f -P(/o)+/(f g{a)<l>(iT-a)da 


D ■ )T 

-^mmy -‘- 


JQ 

cA 


V Pifo)+fo gia)4>ikT-a)da 
f PjM+Ii)' g(a}(t>((i+l)T-a)da 


> 


> 


V Pifo)+fo gia)tt>{iT-a)da 
min{(L)* - Dram) T,Xo + (* + 1) {D* - Dram)T} 

/ P(fo}+f(^ g(a)<f>((i+l)T-a)da 
V P(fo)+f^ g(a)if){kT-a)da 

min{0,Xo + (z + 1) {D* - Dram)T} 

f P{fo)+Io g{a)(l)({i+l)T-a)da 
V PUo)+^Q g{a)4>{kT-a)da 


+ minfc=o,...,i (In 


+ minfc=o,...,i+i In 


which is the hrst inequality (42) with i + 1 in place of i > 0. Furthermore, 0 > {D* — 
Dmm)T + Xi- Combining the previous inequality with definition (39) and the fact that 
Xi+i =Xi + {D* - Dmm)T + Ui, we get 


( P{fo)+!(t g(a)(/)(h+l)r-a)da 


Xi (yD T -\- Ui ^ Vj 

_ maXfc_o,.,.,j+l ^ P(fo)+ff^ g{a)(f>{kT-a)da 

< max{0,xo - (f + 1) (-Dmax - -D*)T} 

, . f]r, f ^Co)+/o^ g{a)<l){{i+l)T-a)da 

k-0,...,i+l g(a)4>{kT-a)da 


Case 3: D* + T 


-1 


which is the second inequality (|42|) with i + 1 in place of i > 0. 

Arguing as in Case 3 in the proof of Claim gives 


Xi 


> Dr 


Xj-)-i Xi (-Djnax D )T -\- Ui 
inequality (|42|, we get 


Combining the previous equality with dehnition 


and 


Xj-|-i Xi {Dmsac D ) T Ui 

< max {0, Xo - i (F>max - D*)T} - (F>max - D*) T 

+ maXfc=o,..., ( In ^ 3iam^^-a)da 


\Pifo)+fo gia)<j>{kT-a)da 


In 


f -P(/o)+/o g{a)4>{{i+l)T-a)da 


V ^{/o)+/o gia)if>{iT-a)da 

< max{- (-Dmax - D*) T,Xo - (i + 1) (-Dmax - D*)T} 

I DM+fo^ g(a)(A(p+l)r-a)d 

’■■■’* V V ^’(/o)+/cf g{a)<j){kT-a)da 

< max{0,Xo - (f + 1) (F>max - A)*)T} 

, f Pifo)+I^ g(aM(i+l)T-a)da'' 

+ max,=0,....+l (^In (, P(/o)+/„%(a)^(.T-a)da . 


which is the second inequality (42) with i + 1 in place of i > 0. Furthermore, Xi — (-Dmax — 


D*)T > 0. Combining the previous inequality with dehnition (39) and the fact that Xj+i = 
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Xi - (-Dmax - D*)T + Ui, we get 


Xi+l Xi (-Djnax D ^ T -\- Ui 

f -P(/o)+/(f gia)(t>iii+l)T-a)da 


> Ui > minfc=o,.... 


,i+l 


V ^’(/o)+/o g{a)<f>{kT-a)da 


> min{0,a;o + (i + 1) (-D* --Dmin)T} 

, f Pifo)+f^ g(a),j)({i+l)T-a)da 

+ niinfc=o,...,i+i I in I —„ ..pA , — r ,— 

V V ^(fo)+Jg g(a)<f>(kT-a)da 

which is the hrst inequality (42) with i + 1 in place of z > 0. This proves Claim 


A.4 Proof of Claim 


Our expressions (34) for the output give 


:.(*) = - (A - A)(« - ^T) + 1„ 


(A.6) 


for all t > 0, where i = [t/T]. We again consider three cases. 


Case 1: Dmin < D* +T ^Xi < U ma x- In this case, (38) and (39) imply that Di = D* + T 


Hence, (A.5)gives 


x{t) = iT)T-^) Xi + In { 

^ ^ ^ ^ ^ \P{fo)+fo g{a)(j)(iT-a)da) 

for all f > 0. The above equality in conjunction with the fact that 0 < 1 — (t — iT)T~^ < 1 


and inequality (42) gives estimates (43)-(44). 

Case 2: D m\r, > D* + T~^Xi. Now our definitions (38) and (39) give Di = Zimin- The 


inequality Zimin > D* +T ^Xi implies that 0 > (Zl* — Zlmin)T +Xi and 0 > Xi, which combined 


with (A.5) give: 


Xi 


Z Pifo)+f^ g(ffl)(^(t-a)da \ ^ < \n ( -^^Zo)+/o'^ g{a)4,{t-a)da 

\P{fo)+f^ g{a)<f>{iT-a)daJ — ^ — V-P(/o)+/cf g{a)(j){iT-a)da 


The above inequality in conjunction with inequalities (42) gives (43)-(44). 

Case 3: D* +T~^Xi > Zlmax- Definitions (38)-(39) imply that Di = Zlmax- The inequality 
D* + T~^Xi > Zlmax implies that x* — (Zlmax — D*)T > 0 and 0 < Xj, which combined with 


(A.5) gives 


1 ( PUo)+Ia g{a)(l>(t-a)da \ , , ^ , -L In Pifo)+In g(a)<^(t-a)da \ 

\P{k)+!^g{a)^{^-a)da) - \P{h)+g(a)^{^-a)da) ' 

which we can combine with (42) to get (43)-(44). This proves Claim 


A.5 Proof of Claim |4| 


Since |ln(x)| 


In (max{x, x ^}) for all x > 0, we get 


In 
In ( 


f P(h)+Jn 9{a')4>{t-aycL \ _ 

V-P(/o)+/o'^ g(a)0(iT-a)da/ 
'max [ g{a)^{t-a)da 

, I -P(/o)+/j^ g{a)cj)(iT-a)da 


P{h)+So g(a)(/)(iT-a)da 1 \ 
Pifo)+f(^ g[a)rf){t-a)da j J 
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and 


= In (max 


f -P(/o)+/o'^ g{a)(j){{i+l)T-a)da P{fo)+f^ g(a)(l>{iT-a)da 
I Pifo)+f(^ g{a)<f>{iT-a)da ’ P{fo)+f^^ g{a)ff>{{i+l)T-a)da 


On the other hand, we can nse (47) to get the following for all i > j: 

< exp(5), 


(A.6) 


(A.7) 


where h = K* ||/o||i exp {—e{iT — A)). Using our bound (46) on 0 and the fact j > [A/T] +1 
(which implies that jT > A, i.e., iT — a > 0 for all a G [0, A]), we get 

P{fo)-K*\\fo\\iexp{-e{{i+l)T-A)) g{a)da ^ P{fo)+f(^ g{a)ff>{{i+l)T-a)da 


< 


-P(/o)+i^*||/o|ll exp(-e(iT-A)) Jg g(a)da 
-P(/o)+-f^*ll/olli exp(-£(p+l)r-A)) g{a)da 
Pifo)-K*\\fo\\iexp(-e{iT-A)) g{a)da 


Pifo)+Io gia)4>iiT-a)da 


and 


P{fo)-K*\\fo\\l exp{-e{t-A)) fg g{a)da 
P{fo)+K*\\fo\\iexp{-6{iT-A)) g{a)da 


< 


■P(/o)+/o gia)4>it-a)da ^ P(/o)+-fC*||/o|li exp(-£(t-A))g(a)da 


Pifo)+f(^ gia)(t>(iT-a)da P{fo)-K*\\fo\\j^exp(-e{iT-A)) gia)da 


for all t > iT. 


Our formulas (35) and (37) for f3 and g imply that /„ g{a)da = 1. Hence, the preceding 


inequalities give 


and 


P{fo)-h ^ P{h)+!o g(a)(/>((i+l)r-a)da ^ P{fo)+h 


< 


P{fo)+h — P{fo)+fg g(a)<l>{iT-a)da — P{fo)-h 


P{fo)-h ^ -P(/o)+/o g(0<^(t-a)da ^ P{fo)+h r ,, ^ > VT 
P{fo)+h — P(/o)+/(f g{a)<l,(iT-a)da — P{fo)-h — ’ 


(A.8) 


where h = K* ||/o|| exp (—^(zT — A)). Combining (A.6) and (A.8), we get: 


\ui\ < In and 


In 


/ P(/o)+/o gia)<l>it-a)da 
V P{fo)+f^ g{a)4>(iT-a)da 


)|<ln(™a) foralH>ff, 


(A.9) 


Using (A.7) and (A.9) we obtain the desired inequality \ui\ < 5 for all i > j. Also, our 
assumptions on k and /o ensure that P{fo) > 0. 


Next, using (47) and (A.9), writing 

P(fo)+h _ 1 I 2h 

P(fo)-h - h-P{fo) ’ 

and using the inequality ln(l + x) < x for all a: > 0, we obtain: 

2h , 2h /;,(exp(5) + l) 


(A.10) 


hi] < 


P(/o) - ft - F(/„) - iglUP/o) ^(/«) 


and 


/ P(/o)+/o g(a)4>{t-a)da 
\P{fo)+f,^ g{a)<f){iT-a)da 


< for all t > iT. 


(A.ll) 
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Since h = K* ||/o||i exp [—e{iT — A)), we obtain 

< ^ ll/o|li(e^W+bexpM) gxp(—eiT) for all i > j and 

^ ^*ll/olli(exp(<5)+l)exp(£^) f 

— Pifo) 1 t. tl ) 

for all i > j and t > iT. This completes the proof of Claim 


f Pifo)+I^ g(a)0(i-a)da 

\P{fo)+f^ g{a)4>{iT-a)da 


(A.12) 

(A.13) 
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